Design of one-dimensional Lambertian diffusers of light 
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Abstract 

We describe a method for designing a one-dimensional random surface that acts 
as a Lambertian diffuser. The method is tested by means of rigorous computer 
simulations and is shown to yield the desired scattering pattern. 
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Optical devices that give rise to a scattered intensity that is proportional 
to the cosine of the scattering angle are frequently used in the optical industry, 
e.g. for calibrating scatterometers [1]. Such diffusers have the property that 
their radiance or luminance is the same in all scattering directions. Due to this 
angular dependence such devices are often referred to as Lambertian diffusers. 
In the visible region of the optical spectrum volume disordered media, e.g. com- 
pacted powdered barium sulphate, and freshly smoked magnesium oxide [2] are 
used as Lambertian diffusers. However, this type of diffuser is inapplicable in 
the infrared region due to its strong absorption and the presence of a specular 
component in the scattered light, in this frequency range. 

The design of a random surface that acts as a Lambertian diffuser, espe- 
cially in the infrared region of the optical spectrum, is therefore a desirable 
goal, and one that has been regarded as difficult to achieve [3]. In this paper 
we present a solution to this problem that is based on an approach used in 
several recent papers to design one-dimensional random surfaces with specified 
scattering properties [4-6], and to fabricate them in the laboratory [5,7]. The 
design of a two-dimensional random surface that acts as a Lambertian diffuser 
will be described elsewhere [8]. 

To motivate the calculations that follow we begin by considering the scatter- 
ing of s-polarized light of frequency uj from a one-dimensional, randomly rough, 
perfectly conducting surface defined by x% = £(xi). The region x% > (,{x\) is 
vacuum, the region x^ < C( x i) is the perfect conductor (Fig. 1). The plane of 
incidence is the xiX3-planc. The surface profile function ((x\) is assumed to be 
a single-valued function of x\ that is differentiable, and to constitute a random 
process. 

The mean differential reflection coefficient (dR/d9 s ), where the angle brack- 
ets denote an average over the ensemble of realizations of the surface profile 
function, is defined such that (dR/d9 s )dd s is the fraction of the total time- 
averaged flux incident on the surface that is scattered into the angular interval 
(9 S ,9 S + d9 s ) in the limit as d9 s — > 0. In the geometrical optics limit of the 
Kirchhoff approximation it is given by [5] 
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In this expression L\ is the length of the Xi-axis covered by the random sur- 
face, 0o an d 0s are the angles of incidence and scattering, respectively, a = 
(u> I c) (cos 0o + cos 9 S ) , and q = (co/c) sin 9 S , k = (co/c) sin 0o . In what follows, we 
will restrict ourselves to the case of normal incidence (0o = 0°), in which case 
Eq. (1) simplifies to 
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Figure 1: The scattering geometry assumed in this paper. 



where a is now given by a = (w/c)(l + cos8 s ). 

We wish to find a surface profile function C(xi) for which the mean differen- 
tial reflection coefficient has the form 
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To this end we write CO^i) m the form [5] 
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Here the {eg} are independent, positive, random deviates, b is a characteristic 
length, and the function s(x\) is defined by [5] 
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where m is a positive integer. Such trapezoidal grooves can be generated exper- 
imentally [5,7]. 

Since the {ce} are positive random deviates, their probability density func- 
tion (pdf) / (7) = (5(7 — a)) is nonzero only for positive values of 7. 

It has been shown [5] that when the surface profile function is given by Eqs. 
(4) and (5), the expression (2) for the mean differential reflection coefficient 
becomes 
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Thus, we find that in the geometrical optics limit of the Kirchhoff approximation 
the mean differential reflection coefficient is determined by the pdf of the coeffi- 
cients {eg} entering the expansion (4), and is independent of the wavelength of 
the incident light. If we make the change of variable tan(# s /2) = jh, < 7/1 < 1, 
so that ±cos6> s = |(1 - j 2 h 2 )/(l + j 2 h 2 ), on combining Eqs. (3) and (6) we 
find that the equation determining f(j) is 

/(-7) + /(7) = 2h (1 1 "ff )2 . (7) 

It follows that 

The preceding results were obtained in the geometrical optics limit of the 
Kirchhoff approximation for a perfectly conducting surface. However, our earlier 
experience in designing surfaces with specified scattering properties [4-6] shows 
that when a surface designed on the basis of these assumptions is ruled on a lossy 
metal, the results of rigorous scattering calculations show that the resulting scat- 
tering pattern retains the form prescribed in the approximate, single-scattering 
calculations. We now demonstrate that such a result is obtained in the context 
of the present problem. 

From the form of f(-f) given in Eq. (8) a long sequence of {q} was generated 
by applying the rejection method [9], and the resulting surface profile function 
C(xi) was generated by the use of Eqs. (4) and (5). We found from numerical 
experiments that in order to have a surface that acts as a Lambertian diffuser 
in reflection the parameter b had to be large. Physically this means that the 
grooves CO^i) have to be wide. 

In Fig. 2 we present the results of rigorous numerical Monte Carlo simulations 
[10] for the angular dependence of the mean differential reflection coefficient 
(dR/d9 s ) for s-polarized incident light of wavelength A = 612. 7nm scattered 
from a randomly rough silver surface of the type described above (noisy curve). 
The value of the dielectric constant of silver at this wavelength is e(w) = —17.2 + 
i0.5. The surface was characterized by the parameters b = 80A = 49^m, h = 0.2, 
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Figure 2: The noisy curve is {dR/d0 s ) calculated by a numerical simulation 
approach for a random silver surface defined by Eqs. (4) and (5) with b = 80A, 
h = 0.2, rn = 1, and the pdf (8), when s-polarized light of wavelength A = 
612. 7nm (e(u>) — — 17.2 -MO. 5) is incident normally on it. The upper solid curve 
is (dR/dd s ) given by Eq. (3). The lower solid curve is the error in the calulated 
mean differential reflection coefficient as measured by its standard deviation. 

and m = 1, and its length used in the simulation was L\ = 164A = 100/im. 
Furthermore, the plot in Fig. 2 was obtained by averaging the results for = 
35, 000 realizations of the surface profile function Ci x i)- Such a large number of 
surface realizations was needed in order to reduce the noise level sufficiently. The 
reason for the slow convergence of the mean DRC with increasing we believe 
is due to the large value of b used in the simulations. Without compromising the 
spatial discretization used in the numerical calculation (Ax\ = 0.164A) needed 
in order to resolve the oscillations of the incident field, only a few grooves s(xi) 
could be included for each realization in the sum (4) defining the surface. 

The lower smooth curve represents an estimate of the error in the calculated 
(dR/d9 s ) due to the use of an finite number of surface realizations for its calcu- 
lation. This error is obtained as the standard deviation of the mean differential 
reflection coefficient (see Ref. [10] for details). 

The upper smooth solid curve in Fig. 2 represents the geometrical optics 
limit of the Kirchhoff approximation, Eq. (3). As can be readily observed from 
this figure, the agreement between the geometrical optics limit of the Kirchhoff 
approximation for a random perfectly conducting surface and the result of rig- 
orous numerical simulations for a real random silver surface is excellent within 
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the noise level. This is indeed the case for all scattering angles S , which we 
find somewhat surprising, since one might have expected the geometrical optics 
approximation to break down for the largest scattering angles. That this is 
not observed in our simulation results is probably an indication that multiple 
scattering processes are of minor importance in the scattering taking place at 
the random surface even for the largest scattering angles. 

Simulations (results not shown) were also performed where the wavelength 
of the incident light was changed by plus and minus 10% from its original value 
of A = 612.7nm. Such changes did not affect the Lambertian nature of the scat- 
tered light in any significant way. This weak wavelength sensitivity is consistent 
with our earlier experience in designing surfaces with specified scattering prop- 
erties [4-6]. Surfaces generated on the basis of different b parameters have also 
been considered. We found that the scattered intensity showed little sensitivity 
to this parameter as long as it is large. 
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